Influence of shear viscosity on the correlation between the triangular flow and initial 

spatial triangularity 



A. K. ChaudhurQ 

Variable Energy Cyclotron Centre, 1/AF, Bidhan Nagar, Kolkata 700 064, India 

(Dated: February 26, 2012) 

In a hydrodynamic model, with fluctuating initial conditions, the correlation between triangular 
flow and initial spatial triangularity is studied. The triangular flow, even in ideal fluid, is only 
weakly correlated with the initial triangularity. The correlation is largely reduced in viscous fluid. 
Elliptic flow on the other hand appears to be strongly correlated with initial eccentricity. Weak 
correlation between triangular flow and initial triangularity indicate that a part of triangular flow is 
unrelated to initial triangularity. Triangularity acquired during the fluid evolution also contributes 
to the triangular flow. 
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In ultra-relativistic nuclear collisions, a dcconfincd 
state of quarks and gluons, commonly called Quark- 
Gluon-Plasma (QGP) is expected to be produced. One 
of the experimental observables of QGP is the azimuthal 
distribution of the produced particles. In a non-zero im- 
pact parameter collision between two identical nuclei, the 
collision zone is asymmetric. Multiple collisions trans- 
form the initial asymmetry into momentum anisotropy. 
Momentum anisotropy is best studied by decomposing it 
in a Fourier series, 
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4> is the azimuthal angle of the detected particle and ip n 
is the plane of the symmetry of initial collision zone. 
For smooth initial matter distribution, plane of sym- 
metry of the collision zone coincides with the reaction 
plane (the plane containing the impact parameter and 
the beam axis), i\) n = ^_Rp,Vn. The odd Fourier co- 
efficients are zero by symmetry. However, fluctuations 
in the positions of the participating nucleons can lead 
to non-smooth density distribution, which will fluctuate 
on event-by-event basis. The participating nucleons then 
determine the symmetry plane (tppp), which fluctuate 
around the reaction plane As a result odd harmon- 
ics, which were exactly zero for smoothed initial distri- 
bution, can be developed. It has been conjectured that 
third hadronic U3, which is response of the initial tri- 
angularity of the medium, is responsible for the observed 
structures in two particle correlation in Au+Au collisions 
0,0, 0|, [|, 0,0. The ridge structure in pp collisions 
also has a natural explanation if odd harmonic flow de- 
velops. Recently, ALICE collaboration has observed odd 
harmonic flows in Pb+Pb collisions |8|. In most central 
collisions, the elliptic flow (1)2) and triangular flow (1)3) 
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are of similar magnitude. In peripheral collisions how- 
ever, elliptic flow dominates. 

In a hydrodynamic model, collective flow is a response 
of the spatial asymmetry of the initial state. For exam- 
ple, elliptic flow is the response of ellipticity of the initial 
medium. If ellipticity in the initial medium is charac- 
terized by spatial eccentricity, £2 = <<p^f2>>' more 
eccentric is the initial medium, more flow will be gener- 
ated, V2 oc e 2 - Similar correlation is expected between 
the triangular flow and initial triangularity. The corre- 
lation between triangular flow and initial triangularity 
is not well studied for fluctuating initial conditions, in 
particular in viscous fluid. Qualitatively, one can argue 
that correlation between momentum anisotropy of pro- 
duced particles and asymmetry of initial density distribu- 
tion will reduce in presence of viscosity. The reasoning is 
simple, in ideal hydrodynamics only a single length scale 
exits in the problem, that of the system size R. Viscos- 
ity introduces an additional length scale, the microscopic 
scattering length. The additional length scale will re- 
duce the correlation between flow coefficients and initial 
spatial asymmetry. In the present brief report, in a hy- 
drodynamic model with fluctuating initial conditions, we 
have studied the effect of (shear) viscosity on the corre- 
lation between triangular flow and initial triangularity. 
For comparison, we have also studied the correlation be- 
tween elliptic flow and initial eccentricity. In ideal fluid, 
elliptic flow shows strong correlation with spatial eccen- 
tricity. The correlation is weakened with viscosity. For 
triangular flow, even in ideal fluid, the correlation is not 
strong. In viscous fluid, it gets even weaker. 

For the fluctuating initial conditions, we have used a 
model of hot spots in the initial states. Similar models 
are used to study elliptic flow in pp collisions at LHC 0- 
[TlT | . In [l2j , a similar model was used to study anisotropy 
in heavy ion collisions due to fluctuating initial condi- 
tions. Recently, we have used the model to study viscous 
effects on elliptic and triangular flow [l3[ . In the model, it 
is assumed that in an impact parameter b collision, each 
participating nucleon pair randomly deposit some energy 
in the reaction volume, and produces a hot spot. The hot 
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spots are assumed to be Gaussian distributed. The initial 
energy density is then super position of N = N participant 
hot spots. 
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is calculated in a 



Glauber model. We also restrict the centre of hotspots 
(r"j) within the transverse area defined by the Glauber 
model of participant distribution. The central density 
£o and the width a are parameters of the model. We 
fix cr=l fm. The central density eo is fixed to reproduce 
approximately the experimental charged particles in a 
peripheral (30-40%) Pb+Pb collisions. 

We characterize the initial density distribution in 
terms of eccentricity C2 and triangularity 63. 
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FIG. 1: In four panels (a)-(d), for fluid viscosity to entropy 
ratio r//s=0, 0.08, 0.12 and 0.16, the correlation between el- 
liptic flow (W2) and initial eccentricity (£2) is shown. Event 
size is N event =500. 
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tp2 and ip3 in Eq l3al3bl are participant plane angle for 
elliptic and triangular flow respectively. Note that for the 
triangularity, we have used the definition due to Teaney 
and Yan Q- An alternate definition was used by Alver 
and Rolland [BJ, where r 3 terms in Eq |3bl are replaced by 
r 2 . However, Teaney and Yan argued from theoretical 
consideration that in the definition of triangularity, r 
terms are more appropriate that r 2 . 

Space-time evolution of the fluid was obtained by solv- 
ing Israel-Stewart's 2nd order theory. We assume that 
in -^5^^=2.76 TeV, Pb+Pb collisions at LHC, a baryon 
free fluid is formed. Only dissipative effect we consider is 
the shear viscosity. Heat conduction and bulk viscosity is 
neglected. Space-time evolution of the fluid was obtained 
by solving the following equations, 
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Eq|4al is the conservation equation for the energy- 
momentum tensor, T^ u = (e + p)u^u v — pg^ w + tt^ , e, p 
and u being the energy density, pressure and fluid veloc- 
ity respectively, ir^ is the shear stress tensor. Eq|4bl 
is the relaxation equation for the shear stress tensor 
n^ v . In Eq |4bl D = u^d a is the convective time deriva- 
tive, V<^u v> = i(Vu" + Vu^) - \{d.u){g^ - u^u v ) 
is a symmetric traceless tensor. 77 is the shear viscos- 
ity and TV is the relaxation time. It may be men- 
tioned that in a conformally symmetric fluid the relax- 
ation equation can contain additional terms 14]. As- 
suming boost-invariance, the equations are solved in 



i In |±* ) coordinates, with the 
code " 'AZHYDRO-KOLKATA" \ developed at the Cy- 
clotron Centre, Kolkata. Details of the code can be found 
in [3. 

Hydrodynamic equations are closed with an equation 
of state (EoS) p — p(e). Currently, there is consensus 
that the confinement-deconfinement transition is a cross 
over and the cross over or the pseudo critical temperature 
for the transition is T c «170 MeV p~6l4T9j . In the present 
study, we use an equation of state where the Wuppertal- 
Budapest [III [l8[ lattice simulations for the deconfined 
phase is smoothly joined at T = T c = 17 '4 MeV, with 
hadronic resonance gas EoS comprising all the resonances 
below mass m res =2.5 GeV. Details of the EoS can be 
found in 

In addition to the initial energy density for which we 
use the model of hot spots, solution of partial differential 
equations (Eqs l4al4bl) requires to specify the fluid veloc- 
ity {v x (x,y),v y (x 7 y)) and shear stress tensor (tt^(x, y)) 
at the initial time Tj. One also need to specify the vis- 
cosity (77) and the relaxation time (tv). A freeze-out pre- 
scription is also needed to convert the information about 
fluid energy density and velocity to particle spectra. We 
assume that the fluid is thermalized at r^=0.6 fm and the 
initial fluid velocity is zero, v x (x,y) = v y (x,y) = 0. We 
initialize the shear stress tensor to boost-invariant val- 
ues, tt xx — tt vv — 2ri/3Ti, n xy —0 and for the relaxation 
time, we use the Boltzmann estimate tv = 3rj/2p. We 
also assume that the viscosity to entropy density (77/ s) 
remains a constant throughout the evolution and simu- 
late Pb+Pb collisions for a range of r//s. The freeze-out 
is fixed at T F =130 MeV. 

For fluid viscosity to entropy ratio rj/s=0, 0.08, 0.12 
and 0.16, we have simulated b=8.9 fm Pb+Pb collisions. 
b=8.9 fm collisions approximately corresponds to 30-40% 
collision. In viscous evolution, entropy is generated. To 
account for the entropy generation, the Gaussian den- 
sity £0 was reduced with increasing viscosity, such that 
in ideal and viscous fluid, on the average, 7r~ multiplic- 
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FIG. 2: same as in Fig[T] but for the triangular flow(«3) and 
initial triangularity (£3). 



ity remains the same. In the present study, we have used 
N event =500 events. In each event, Israel-Stewart's hy- 
drodynamic equations are solved and from the freeze-out 
surface, invariant distribution ( dy<$p T ) ^ or 7r_ was 
tained. In analogy to Eq |3al3bl invariant distribution 
can be characterized by 'harmonic flow coefficients' [2lj |. 



f P Tdp T d(j)e m ' t > ,, dN . 



In a boost-invariant version of hydrodynamics, flow co- 
efficients are rapidity independent and in the following, 
we drop the rapidity dependence. Present simulations 
are applicable only in the central rapidity region, y «0, 
where boost-invariance is most justified. 



TABLE I: Correlation measure for elliptic and triangular flow, 
as a function of viscosity over entropy ratio. 





?=0 


2=0.08 


2=0.12 


2=0.16 


C measure. (^2 ) 


0.46 
8.77 


0.53 
s.77 


0.92 
S.77 


1.77 
s.77 


Cmeasure (^3) 


2.46 
R.77 


2.95 
s.77 


3.91 
S.77 


4.50 
S.77 





In Fig|TJ in 4 panels, for fluid viscosity 77/5=0, 0.08, 
0.12 and 0.16, simulated elliptic flow (7J2) is plotted 
against initial eccentricity (£2). Each panel contains 500 
data points. In ideal fluid, v-i and £2 are strongly cor- 
related, «2 oc £2. Evidently, the correlation is gradually 
weakened as the viscosity of the fluid is increased. The 
result is not unexpected. As argued earlier, in viscous 
fluid, correlation between elliptic flow and initial eccen- 
tricity is reduced due to introduction of the additional 
length scale. Correlation between triangular flow (7J3) 
and initial triangularity (£3) is studied in FigJ^l The re- 
sults are more interesting. Even in ideal fluid, v$ and £3 
are not strongly correlated. In viscous fluid correlations 
are even worse. Indeed, for fluid viscosity ri/s=0.12-0.16, 
it appears that flow coefficients are marginally related to 
initial triangularity. 



Qualitatively, from FigQ] and [2 one understand that 
the correlation between triangular flow and initial tri- 
angularity is much less than the correlation between the 
elliptic flow and initial eccentricity. One also understands 
that the correlation reduces with viscosity. We can ob- 
tain a quantitative measure of the correlation between 
flow coefficients (W2/3) an d initial asymmetry parameter 
(£2/3). We note that for a perfect correlation, v n oc e n 
and simulated flow coefficients will fall on a straight line. 
Dispersion of the flow coefficients around the best fitted 
straight line then gives a measure of the correlation. We 
thus define a correlation measure function C measure , 
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which measure the dispersion of the simulated flow coef- 
ficients from a best fitted strai ght line, relative to com- 
pletely random flow coefficients |22[ . In order to compare 
the correlation between V2 and £2 and the correlation 
between i> 3 and £3, we also rescale the flow coefficients 
(^2/3) and asymmetry parameters (£2/3) to vary between 
(0-1), such that the dispersion is measured in a com- 
mon SCale. Cmeasure 

varies between 0-1. If simulated 
flows are perfectly correlated with the asymmetry mea- 
sure, C measure = 0, at the opposite limit, when they are 
perfectly uncorrelated (random) C measure = 1. In ta- 
bleU we have noted Cmeasure for elliptic and triangular 
flow, as a function of 77/s. With increasing fluid viscos- 
ity, correlation between the flow coefficient and spatial 
asymmetry parameter is reduced. In ideal and mini- 
mally viscous fluid, elliptic flow and initial eccentricity 
are mostly correlated, C m easure{v2) ~ 0.05. For more 
viscous fluid 77/s=0.12-0.16, correlation though reduced, 
remains strong, C mea sure(v2) ~ 0.1 — 0.2. If we inter- 
pret C m easure as the fraction of flow unrelated to the 
initial spatial asymmetry, for elliptic flow, the fraction 
is small, less than ~10% for fluid viscosity over entropy 
ratio ?7/s=0-0.12. Triangular flow on the other hand ap- 
pears to be highly uncorrelated. Even for ideal fluid, 
Cmeasureivs) pa 0.3. In the above interpretation, ~30% 
of triangular flow is unrelated to the initial triangularity. 
For more viscous fluid, ?7/s=0.12-0.16, the fraction is in- 
creased to ~50%. Comparatively large value of C measure 
for triangular flow raises an important question. Is tri- 
angular flow is response of the initial triangularity of the 
medium only? One may argue that in RHIC/LHC en- 
ergy collisions, viscosity of the produced fluid is not large 
and the simulation results for fluid viscosity ?y/s=0.12- 
0.16 will not be of any practical concern. However the 
correlation between V3 and £3 is still weak in ideal or 
minimally viscous fluid. While for elliptic flow, ^95% of 
the flow is related to initial spatial asymmetry, for the 
triangular flow, the fraction is only ~65-70%. It is rea- 
sonable to conjecture that a large part of the triangular 
flow is unrelated to the initial triangularity. What are the 
mechanisms by which the system acquires triangularity 
is uncertain. Processes like jet quenching, Cerenkov ra- 
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diation etc. may introduce triangularity in the system. 
Indeed, hydrodynamical simulations of jet quenching do 
indicate development of triangularity in density distribu- 
tion at late time (see Fig. 2 of ref . [23j ) . 

To summarize, in a hydrodynamic model, with fluctu- 
ating initial conditions, we have studied the correlation 
between triangular flow and initial triangularity of the 
medium. In ideal or minimally viscous fluid, triangular 
flow is only weakly correlated with initial triangularity. 



In more viscous fluid, the correlation gets even weaker. 
Elliptic flow on the other hand is strongly correlated with 
initial eccentricity in ideal or viscous fluid. Weak corre- 
lation between triangular flow and initial triangularity 
strongly indicate that a part of the triangular flow is 
unrelated to initial triangularity of the medium. Final 
state triangularity, generated by unknown mechanisms, 
also contributes to the triangular flow. 
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